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Abstract 

We derive a mathematical model for eddy currents in two dimensional geometries 
where the conductors are thin domains. We assume that the current flows in the xz- 
direction and the inductors are domains with small diameters of order 0(e). The model 
is derived by taking the limit e — > 0. A convergence rate of 0(e a ) with < a < 1/2 in 
the L 2 -norm is shown as well as weak convergence in the W x ' v spaces for 1 < p < 2. 

1 Introduction 

Mathematical modelling of eddy current problems often involves multiple conductors with 
various sizes. Typically, electrotechnical devices involve thin conductors as wires or coils as 
well as massive conductors. Numerical solution of such problems may then encounter serious 
difficulties in the choice of the domain meshes which can in particular lead to ill conditioning. 
Asymptotic analysis of these problems appears as an efficient tool to obtain limit problems 
that are simpler to solve and better conditioned. 

We consider, in the present work, a two-dimensional eddy current problem, formulated in 
terms of a scalar potential in the whole plane. The electrically conducting domain consists 
in a "thick" conductor and two "thin" domains assumed to carry the same current with 
opposite sign. In terms of the current conservation principle, this means that these inductors 
are assumed to be virtually linked at the infinity. The derivation of the model for thin 
inductors is obtained by assuming that these domains are of small diameters of order e <C 1 . 
We show that taking the limit when these diameters tend to zero leads to a singular elliptic 
problem, the singularity being due to the presence of Dirac measures. 

The outline of the paper is the following: We start in Section 2 by deriving the considered eddy 
current model from a 3-D model. We emphasize on a careful modelling that takes into account 
the total current flowing in the inductors. In Section 3, we state the main convergence result 
and prove it through some preliminary lemmas. Section 4 is devoted to further convergence 
results in W 1,p spaces. 

2 Statement of the problem 

Let A = ft x M denote a cylindrical conductor where f2 is a domain in K 2 with a smooth 
boundary T. We assume that the domain Q is the union of three connected domains flk with 



respective boundaries Tk, k = 0, 1, 2 (see Figure 1), and that the closures of the domains f2fc 
are disjointed. We shall also deal with the complement O' = R 2 \ f2 of f2. 
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Figure 1: A typical configuration of the conductors 

In the following we shall make use of generic constants that do not depend on the small 
parameter e. Time harmonic eddy currents equations read: 



curl H J = 
curl H = 

iojfxH + curl (cr _1 J) 
div (fiH) = 



in A, 

in R 3 \ A, 
in A, 
in R 3 . 



(2.1) 



Here the vector fields H and J denote respectively the magnetic field and the current density. 
Moreover, a and /x are respectively the electric conductivity and the magnetic permeability. 
We assume for the sake of simplicity that a and fi are positive constants. In order to take 
advantage of the geometry of A, we seek unknowns in the form: 

H(x 1 ,x 2 ,x 3 ) = H 1 (x 1 ,x 2 )e 1 + H 2 {x 1 , x 2 ), 0) e 2 , 
J(xi,x 2 ,xz) = J(xi,x 2 )e 3 , 

where (ei, e 2 , e 3 ) is the canonical basis of R 3 . Equations (|2.1() become then 



curllf-J = in fi, 

curlJf = in Q' , 

iuj/iH + curl (<7 -1 J) = in f2, 

div H = in R 2 . 



(2.2) 
(2.3) 
(2.4) 
(2.5) 



where curl and curl denote respectively the scalar and vector curl operator in 2-D, i.e. 

dH 2 dHi 



curl u 



dxi dx 2 



dip dip 
curl ip = - — ei - - — e 2 . 

ox 2 ax\ 



Our aim now is to derive a simple model for eddy currents. Using Equation l|2.5l) . we deduce 
the existence of a scalar potential u : R 2 — > C such that 



(xH = curl u 



(2.6) 



Equation l|2.2|l yields 



curl curl u 
curl curl u 



fi J 




in Q, 

in n', 



2 



or equivalently, 

—Am = iij in Q, 

2.7 

Au = in 0'. 

On the other hand we obtain from l|2.4|l and H2.6|l . 

curl (iivu + a~ l J) = 0. 

Whence 

iujou + J = aCk in (2-8) 
where Ct are complex constants, for k = 0, 1, 2. Replacing this in i|2.7[) . we obtain 

— Au + iujjj.au = jiaCk in f2&, fc = 0, 1, 2, 
Ait = in Q.'. 

Finally, various considerations dealing with interface conditions and the behaviour at the 
infinity lead to the problem: 

— Au + iujjj.au = jjaCk in k = 0, 1, 2, 
Ait = in fi', 

N = on r, , 2 g ^ 



9u 



on r, 

as \x\ — » oo. 



9n 

u(a;) =a + 0(|x|- 1 ) 



Here above, [ • ] denotes the jump of a function across the boundary T, this jump being equal 
to the external trace minus the internal one. It remains to determine the constants Ck in 
function of problem data. For this end, it turns out to be realistic to prescribe the total 
current in each conductor. We then assume that this quantity, denoted by I is given as 

Jdx = - J dx = I. (2.10) 

Sli JQ.2 

Note that the first identity is imposed in order to enforce a current conservation principle. 
For the same reason, we impose 

Jdx = 0. 



Making use of these conditions, we obtain for the constants Cfc, the values 

Co = iuJUx, 

„ I 

C% = iuju\ + 



C 2 = IUJU2 



a\^\ 
I 



where \Qk\ stands for the measure of f2fc and Uk is the average of u on Qk, i-e. 

™k ■= 77TT / udx. 
"* in. 
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We obtain the problem: 



— Au + iLu/ia(u — uq) = in tto, 

— Au + iixi[ia{u — u\) = j^r- r in Qi, 

l"i| 

— Au + iu!fj,a(u - u 2 ) = -7777 in ^2, 

Am = in f2', 

[u] = on r, 

du 

dn 

u{x) ^a + OUx]- 1 ) 



(2.11) 



= on r, 

as \x\ — > 00. 



Note that, owing to (|2.1U|) . the solution of Problem (|2.11|) is known up to an additive constant. 
For this reason, we impose the condition 

u = 0, (2.12) 

which enforces a value for the constant a. 

Let us prove that Problem 1)2. 11[) - H2. 1211 has a unique solution. We define, for this end, the 
Beppo-Levi space (see [5]), 

W 1 ^ 2 ) := ju; pv £ L 2 (R 2 ), Vw £ L 2 (R 2 ) 2 }, 

where p is the weight function given by 

P[X) = (l + M)log(2 +N) - < 2 ' 13 > 
We furthermore define the space 

V :={ve W^M 2 ); vo = 0}. 
It is well known (cf. ^UJ) that the semi- norm 



r 1 
( / \Vv\ 2 dxY 



is a norm on the space V, equivalent to the one induced by W (K ), i.e. we have in particular 

\\pv\\ L 2 (M 2 ) < C \v\ w i (M 2) V v G V. (2-14) 
Here and in the following |Vu| stands for the function 



/ dv dv dv dv \i 
\dxi dxi dx 2 dx 2 ) 

A variational formulation of (|2.11l) consists in seeking a function u 6 V such that 

2 

/ Vu-Vvdx + (u -u k )vdx = filivx -v 2 ) VueV, (2.15) 

where /? = w/xc and v is the complex conjugate of v. 
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Theorem 2.1. Problem H2.15[) has a unique solution. 

Proof. Let us define, for u, v G V, the scsquilincar and antilincar forms, 

2 

a(u, v) := / Vu • Vw dx + i/3 >^ / (u—Uk)vdx, 
:= fil(vi - v 2 ). 

The forms a and L are obviously continuous. In addition, since 

(u — Uk)vdx— / (u — Uk)(v — vt) dx, 
then we have 

r 2 r 

v k )v dx 

k=0' " 
2 

\v — vt\ 2 dx. 



2 

a(v,v)= / |Vu| 2 dx + i/3 V] / (u 

JR 2 fc=Q 7n fc 

2 

= / |Vu| 2 dx + i/3V / 
•/r 2 t._n-/n 



Then 



Re(a(v,v)) = / |Vw| 2 dx. 
Jr 2 



We deduce then that a is coercive on V and the Lax-Milgram theorem gives the existence 
and uniqueness of a solution iiGl^to i|2.15|) . □ 

We now consider that the domains fii and O2 are thin in the following sense: we define the 
domain 0| := f2^ by 

ni = z k + eh k fc = l,2, 

where e is a small positive number, z k £ R 2 , and tt k is a smooth domain in M 2 . We assume 
furthermore that the domains fl k and Slo are disjointed for e small enough. Furthermore, we 
denote in the following by f2 e the union Slo Uf2f Uil|- Finally, let us mention that, throughout 
this paper, C, Ci, C2, . . . will stand for generic constants that do not depend on e. Our aim 
is to study the asymptotic behavior, as e — » 0, of the solution it to Problem l|2.11|l . 

3 The limit problem 

Let us first, for clarity, rewrite Problem (|2.11() with the parameter e. Denoting by xo and x% 
the characteristic functions of f2o and f2i , respectively, we have 



2 e e 

k=l ' 1 ' ' 2 ' 



in R 2 , 



(3.1) 

u e (x) = a + 0(|x| _1 ) as |x| — > 00. 

Let us recall that the condition (|2.12(l fixes the value of a. 

We next define the weighted space that will be used for convergence results: 

L 2 (R 2 ) = {v; pv S L 2 (R 2 )}. 
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We also define a problem that will be defined as the limit problem. This one is the following: 

- Au + iPxou = fJ,I(8 Zl -8 Za ) inR 2 , 

l v J 

u(x) = a + 0(\x\~ ) as \x\ — > oo, 

where S Zk is the Dirac measure concentrated at z^. For Problem (|3.2() . we need a uniqueness 
result. Let us define for this the notion of weak solution. We shall say in the sequel that u is 
a weak L 2 -solution of Problem (|3.2(l if u G L 2 (R 2 ) and if we have 



u(-ATp + if3xo^p)dx = ^I(Tp( Zl )-Tp(z 2 )) V^e^R 2 ), (3.3) 

where £^(R 2 ) is the space of indefinitely differentiable functions with compact support in R 2 . 
Lemma 3.1. Problem (|3.2[) has at most one weak L 2 -solution. 

Proof. Let ui and u 2 denote two weak L 2 -solutions of (|3.3|l . The difference u = u\ — u 2 
satisfies then 

ui-ATp + ifixo^dx = V</?6^(R 2 ). 



This relation is still true for all functions <p € L 2 (R 2 ) with 

p 2 ulpdx = V t/> e L 2 (R 2 ), 



where 

-A7p + if3xo~(p = p 2 Jj inR 2 . (3.4) 
Note that Equation l|3.4() admits a unique solution in T / l /1 (R 2 ). Choosing ip — u, we deduce 



/ p 2 \u\ 2 dx = 0. 
Jr 2 



This implies u = and uniqueness follows. □ 
We now state the first convergence result. 

Theorem 3.1. The sequence (u e ) converges in L 2 (R 2 ), when e — > 0, to the unique solution 
of Problem 

The remaining of this section is devoted to the proof of Theorem 13.11 It is clear that the 
structure of the right-hand side in Problem (|3.1|1 suggests that the convergence cannot be 
obtained in the space T / l /1 (R 2 ). To obtain a weaker result we resort to a duality technique 
due to Lions-Magenes (E], p. 177) and Damlamian-Ta Tsien Li 

Let, in the following, B denote a ball that contains the domains Slo , and fl 2 for all e«l. 
Multiplying Equation (|3.1(l by a test function cp 6 V Ci H 2 oc (R 2 ) , and using the Green formula, 
we obtain 

2 

— / u e Aip dx + i(3 / (u e — u%)Tpdx + i(3 / u e 7p dx = pl(ip 1 — ip 2 )- 
JR 2 k=1 JQ.% Jn a 



Since 



(u e -u%)cpdx = u e ((p-p k )dx= (u £ -u%)(tp - ip 6 k )dx, (3.5) 
J at Jnt 
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where (p e is the solution in V (~l H 2 oc (R 2 ) of 

2 



fe=l 

Lemma 3.2. We have the estimates: 



fc=i 



(3.6) 



we deduce that 

r 2 _ 

/ " e ( - ATp + iP^xUv ~ Vk) + i(3xo^p) dx = pl(ip 1 - <p 2 ). 
Jr2 k=i 

Let ip denote a function in L 2 (R 2 ). Identity .(j|> can also be written as 

p 2 u e ipdx = fil((fil - tfy, (3.7) 



-Aif e +i(3j2xl(^-^ k )+i0Xo^ = p 2 ^ inM 2 . (3.8) 



\\V(p e \\ L 2 {R2) 2 + \\tp € \\ L 2 {no) +e 1 ^2\\<p e - VfclU=(nj) < C \\fyip\\ L 2 (v) . (3.9) 



W\\h*{b) <0\\fnP\\ L 2 m , (3.10) 
for each ball B of R 2 containing Vt e . 

Proof. By the Green's formula, we have from (|3.8|1 and Identity l|3.5|l 

[ \Vip e \ 2 dx + if3Y] [ \ip e - fii\ 2 dx + i/3 [ \ip e \ 2 dx = [ p 2 ^ dx. 
From this and (|2.14|l we deduce that 

/ \WLp E \ 2 dx< 11^-01^2(82)1^^1^2(82) < C||pVIIl2(R2) ||V(^ £ || L 2( R 2)2, 
JR 2 

and then 

I 

\V^\ 2 d x y <C\\(nP\\ L , m . (3.11) 

Therefore, the sequence ((f e ) is bounded in M /1 (R 2 ). The L 2 -error estimate is obtained by 
using the Poincare-Wirtinger inequality (see [Q, p. 194). We have indeed by using H3.11J) . and 
since the diameter of 0| is an 0(e), 

II^II^COo) < Cl ||V^ e || L 2 (R 2)2 < C 2 ||p-0|| L 2( K 2), 

11^ - <Pk\\L 2 (n> k ) < C 3 e ||V^ £ || L 2 (R 2)2 < C 4 e \\p^\\l 2 (r 2 ), k — 1,2. 

In order to prove the iJ 2 -estimate, we use standard regularity results for elliptic equations 
(See 0, P- 183 for instance). We obtain for any ball B of R 2 containing fl e , and any regular 
domain D containing B, 

2 

\\<p e \\ir>(B) < Ci (\\<P e \\m(D) + \\p 2 M\l 2 (d) + y'h 2 ^) +^2\\<P e - &i\\L 2 {si' k )) 

k=l 

< C 2 \\P^\\l 2 (R 2 )- 

Note that the constant C2 depends on the domain B but does not depend on e. □ 
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The estimates obtained in Lemma l3~21 enable concluding that a subsequence of (p e ) converges 
toward tp weakly in H 2 (B) for any ball B of K 2 . We now characterize the limit function. 

Lemma 3.3. The sequence (ip e ) converges, when e — * 0, in W /1 (R 2 ) to the unique solution 
of the equation: 

-Aip + if3 Xo p = p 2 i> inM 2 , (3.12) 
Moreover, we have the error estimates 

2 

X] H^ C ~ ^Hi 2 (^fc) + 11^ ~ ^Wh'^b) + ||V(^ £ - <p)\\l*(b?)' < Ce ||WIIl2(r2), (3.13) 

k=l 

for any ball B of R 2 containing £! e . 

Proof. Let (f) e = ip e — (p. Then (ff e V (1 H^ a C (K 2 ) and satisfies the variational equation 

2 

/ Vcjf ■ Vvdx + i{3 / frvdx + ipy^ (<P e ~ filP dx = V v € V. 
Choosing v = (jf , we obtain 

2 

/ \V(j) e \ 2 dx + if3 [ \ct) e \ 2 dx = -i(3j2 f ~ ffiW dx - 

JW? JQ.0 k= i Jto% 

Then using the estimates (|3.9fl . we have 

[ \V<f> e \ 2 dx+ f W\ 2 dx< 13^(1 \<P e ~ V%\ 2 dx)* ( f W\ 2 dxY 

< Cl£ ||/3V||i3(R2) ||V0 £ || L 2 (K 2)2. 

Therefore, we have the bounds 

||V0 e || i 2 (H2) 2 < Cie ||pV|U=(R=), (3.14) 
ll^ll^(no) < C 2 ||V^|U2 (R2) 2 <C 3 e\\piP\\ L * m . (3.15) 
The sequence (y> £ ) converges then to ip strongly in W /1 (IR 2 ), which yields the limit problem 

irm . 

To prove the L 2 -error estimate, we have from (I2.14[) and l|3.14[l . for k = 1, 2, 

W\\l*{q.d < Ci ||p0 e |U 2 (nj) < Ci ||p</> c || L 2 (R 2) < C 2 ||V<?!) e || L 2( R 2 ) < C 3 e HpVlU 2 ^)- 
The ff 2 -estimate is handled in the following way: By subtracting l|3.8[) from l|3.12|l . we obtain 



2 



-A0 e = -z/3 X o0 e -i/3^Xfc(^-^) in 



k=l 

Using 13.9|l . I|3.15() and classical regularity results for elliptic problems (See P- 183 for 
instance), we get 

2 

W\\h 2 (b) < C x (\\4?\\hhd) + ll#ll^(Oo) + I] llv> e - Kh^n 

fe=i 

< C 2 U ||V0 £ || L 2( R 2 ) 2 + ||0 £ ||L2(O n ) + 6 ||p^||i2 (R 2 ) 2 

< C 3 e||pV||i2(K 2 ), 
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for all compact subsets B of K 2 and all regular domains D that contain B. Note that the 
constant C depends actually on B. □ 

We are now ready to obtain the first convergence result for u e . 

Theorem 3.2. There exists a constant C, independent of e, such that 

||p( W -u £ )|| L 2 (R2) <Ce a / 2 0<a<l, 

Proof. Consider the problem 1)3. 7)1 and the following one, for ip G L 2 (R 2 ), 



p 2 u tpdx = pi (ip(zi) - ip(z 2 ))- 



(3.16) 



where tp is the solution of Problem 1)3.12)1 . Then 



1 



\n 



p (u e - u)vp dx = pi [ —— ip e dx - tp{zi) ) - pi [ —— ip e dx - tp(z 2 ) ) ■ (3.17) 



ll Jn\ 



1 



n 



2 1 Jn 



Since tp G H 2 (B) C C°' a (B) for all a with < a < 1 (see £Q for instance) and all compact 
subsets B of R 2 , we have for k = 1,2, 



fel Jn 



tp(x) dx - p(z k ] 



< 



in 



|^(x) - <^(z fe )| efe 



< Ce a . 



(3.18) 



Furthermore, we have from 1)3. 13)1 . the imbedding H 2 (B) c C°(B) and the mean value 
theorem, 



1 



05 



(<p e -tp)dx <C 1 \\p t - p\\ c o {B) 



< C 2 \\<p e - <p\\h*(b) 

< C 3 e \\pip\\mm)- 



(3.19) 



Recalling (|3~T7| and using IpTTSj) . (|3.19[l . we get 



lim / (u e - u) p 2 ^ = Vi/)eL 2 ,(l 2 ). 

The sequence (w e ) converges then weakly to u in L 2 (R 2 ). To obtain the strong convergence 
of u e , we choose ip = (u e - u) G L 2 (M. 2 ) in ipDTjl . We have by using again (|3~T%|) . {HUSt, 



I|P(W £ - W)|||2 (R 2' 



2 2 ^ 

-^ 7 X] / i^-^dx +plY] it^tj / ^da; 
*=i 1 •'"S fe=i 1 kl Jn i 



<C A e + C 5 e a < Ce a . 



□ 
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4 Sharper convergence results 

The convergence result obtained in the previous section can be improved, as we shall show 
hereafter, by using the technique of renormalized solutions for elliptic equations following 
Boccardo - Gallouet [2] and Murat jZj. To simplify the settings, we shall sometimes resort to 
writing Problem 1)3. 1|) as a system of two coupled equations involving real valued unknowns. 
Let us denote, for a complex number z, by zr and z\ its real and imaginary parts respectively. 
Equation i|3.1[l can be written: 



2 e £ 

k=l ' X ' ' 2 ' 



inK 2 , (4.1) 



-Auf + /3^Xfc(«fl-«fc,B)+/?Xo«k = inM 2 , (4.2) 
fc=i 

w j?( x ) = a R + ^(M _1 ) as M ~~ * °°j (4-3) 

uf (x) = a/ + 0(|x| _1 ) as \x\ — > oo. (4.4) 

We start by deriving L 2 and L 1 uniform estimates. 

Lemma 4.1. We have the estimates: 

2 

||p W e || L 2 (R2) + e-5 ^ - ul\\ L 2 (ni) < C, (4.5) 
fe=i 

2 

|| P V|| il(K2) + e-i H" £ " "kll^Cnj) < C (4.6) 
fe=i 

Proof. The estimate on Hp w e || ^,2 (m 2 ) is obtained from Theorem 13.21 and from the fact that 
pu e L 2 (R 2 ). Next, The Holder's inequality gives 

/ p 2 \u e \dx<([ \pu £ \ 2 dx) ( I p 2 dx) <C 2 ( I \pu e \ 2 dx 
Jr 2 \Jr 2 J \Jw 2 J \Jm 2 

Using a variational formulation of Problem l|3.1|) , we then obtain the bound 

r 2 1 

Jr 2 |S2 fc l 

The Poincare-Wirtinger inequality yields for k — 1,2, 



|u e -ulYdx < Cie A / \Vu £ \ z dx <C 2 e- 
Jm 2 

Again, the Cauchy-Schwarz inequality gives the L 1 -estimate: 



\u e - u%\dx < I \u e -ulfdx)' <Ce^. 



□ 

We now need a technical result before proving a convergence result. The result, which is a 
variant of the Poincare-Wirtinger inequality, can be established by an analogous proof. 
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Lemma 4.2. There exists a constant C such that 

\Mlv{b) < C\\Vv\\ LP{B) 2 V v € W^iB) with / udx = 0, 

Jf2 

where 1 < p < oo and £? is any compact subset of R 2 that contains O c . 

Theorem 4.1. The sequence (u e ) converges weakly in W lp (B), 1 < p < 2, toward the 
unique solution u of Problem 13.2(1 in each ball B containing fl . 

Proof. For an integer m, we define a subset B m of -B by 

= {x E B- 2 m < max{\u R (x)\, \u}(x)\} < 2 m+1 }. 

Let ip m stand for the truncature function defined by 

'0 if 0<s<2 m , 

ipm(s) := { s - 2 m if 2 m < s < 2 m+1 , 
2 m if 2 m+1 < s, 

extended to R by oddity. Multiplying Equation l|4.1[) by ip m (u R ) and Equation (|4.2|l by 
^mtw}), integrating on R 2 , using the Green formula and summing up, we get 

^' m (u R )\Wu R \ 2 dx+ I *P' m (u})\Vu}\ 2 dx-f3 I X ou^ m (u e R ) dx 
it 2 ji 2 

2 



+ P xou^miu}) dx - pY] / Xfe«-Wj,fe)^m(«fl)rfa; 

JR 2 fc=l^ R2 
2 

/ Xk( u R - UR,k)*P™(u e I ) dx = r e iP m (u e R )dx, (4.7) 

t~! JR 2 JR 2 



fe=i 

where 



Note that we have 

||r e ]U I(K 3) < C. (4.8) 
Since ip' m > and |?/> m (u e )| < 2 m , we have by using l|4.6fl. 



i / |V,'|^/.r : C. (4<)) 



2 



where C is independent of e and to. Let p denote a real number with 1 < p < 2. We have 
from the Holder inequality 

J \Vu e \ p dx<(J |Vu £ | 2 dx) "IB'J 1 -^. (4.10) 

Since |u e | > 2 m on B m , we have by using the Holder inequality, 

\BL\<J- f \ u *\dx<±-( [ \u £ \'dxY\m 



om / 1 — On_ , , 
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for all s, s' > 1 with 1/s + 1/s' = 1. Hence 
Using I4.10|l and l|4.H(l yields then 



B jv^rdx< 2m(s( J /2) . p/2) (/ B j,T^ 

We choose here s > p/(2 — p) so that s(l — p/2) — p/2 > 0. Therefore 

E 1 i v "'i Ptfa ^ g E 2m( ,(i-p/ 2 )-p/ 2 ) (/ fl . i«-r^) 1_J - ( 4 -n) 

m>0 m m > m 

From the discrete Holder inequality 

E a »^<(E<T(E<f forry>i, ~ + i = i, 

m mm 

Inequality l|4.11|l yields 

E 1 iVnT ^<g(E 2 ^ (1 -t/ 2 )- P / 2 ) )"(E(/ £ 



m>0 m>0 

Choosing r' = 2/(2 — p), we obtain 



\u- 



r'(l-p/2)x ^ 



E / |Vu e | p ^<c(E / IwTda;) 1 (4.12) 



We next define 

B f = {i6B; < max{|ufj(x)|, |uf(x)|} < 1}, 

which clearly implies 

B = & U ( |J 

m>0 

In order to estimate u e in W 1,p (B e ), we define the truncation function 

!1 if s > 1 
s if - 1 < s < 1 
-1 if s < -1. 

Multiplying Equation Ij4.1jl by T(u R ), Equation (|4.2|) by T(u}), integrating on R 2 , using the 
Green formula and summing up, we obtain 

/ T\u R )\Vu R \ 2 dx+ f T'(u e I )\Vu e I \ 2 dx-P [ X ou}T(u R )dx 
Jm 2 Jr. 2 Jr 2 

2 

+ /?/ X vu R T(u})dx-pJ2 [ xl(u}-u}^T(u R )dx 
Jm 2 k=1 Jr 2 

2 

+ /?E / X%{u R -u R<k )T{u\)dx= f r*T{u R )dx. 
k=1 Jr 2 Jr 2 
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Using I4.6fl . the bound l|4.8f) and the properties \T(s)\ < 1, V > 0, we deduce 
[ \Vu' R \ 2 dx+j \Vu\\ 2 dx< ||r e || Ll(R2) +/3(||^|| L i (0o) + |K|| L i(o ) 

2 2 

+ W U R ~ ^R-khnni) + hi - u} >k \\ L Kn% 

fc=l k=l 

< c. 

This yields 

f \S7u e \ p dx < C. (4.13) 
Jb> 

Combining l|4.13[l and 14.121) . we have then in particular 

J ^u^dx < C ^1+ (y I^T^) 1 fors> ^/ ( 4 - 14 ) 
We use successively the Gagliardo-Nirenberg (see Friedman |3J, p. 27) and Lemma f4. 21 to get 

u e \ s dxY <C( [ \Vu e \ p dxy([ \u e \dx Xl 

B J B J B 

with < A < 1 and such that 

1-i 



A - 3 i 
2 p 



Using l|4.6[> yields 



B 



where C = C(B). Whence, from (|4.14l) . 

A 3 (l-p/2) ' 



J \Vu e \ p dx <C (l+ (J \Vu e \ p dx^j " J, (4.15) 



for all s > p/(2 — p) and < A < 1. Let us choose for s the value (1 +p)/(2 — p) that yields 

As ( 1 _P\ _ 2 ~P 
p\ 2/ 3p-2 

We then deduce from Ij4.15|l the bound 



■( 1 -f)-S <1 fo rl<p<2. 



|Vu £ | p rfa; < C, 

B 

with C — C(B). Therefore, the sequence (u e ) is bounded in W 1,P (B) for all balls B that 
contain Q . From this, we deduce that a subsequence of (u e ), still denoted by (u e ), satisfies 

u e ^-u* in W 1 ' P (B). 

From the compactness of the imbedding W 1,P (B) C L q (B) for 1 < g < 2p/(2 —p), we have 

u e -> u* in L«(B) for 1 < q < 

2-p 
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Theorem 13.21 implies u* — u. Thus, the subsequence of (u e ) converges strongly to u in L q (B). 
Let us show that the convergence to the solution of (|3.2H takes place in W 1,p (B)-weak for all 
bounded balls B of R 2 . We have from (|3.1() for all <p € W 1,p (B) extended by zero outside B, 
with l/p + l/p' = 1, 

We have 

Vu £ -VTpdx -> I Vu-VTpdx. (4.16) 



Next, using l|4.5[) . we have for k = 1,2, 

(w e - u%)Xpdx < \\u € - ul\\ L 2 {ni) \\cp\\ L 2 {Qi) < Ce? |M|z,2 (n e 



Therefore 

/ {u e -u%)lpdx ->0. (4.17) 
For the term involving Oo, we deduce from Lemma 14.11 

u e 7pdx—> / uipdx. (4-18) 



Finally, since p 1 > 2, then we have the imbedding of W l p (B) into C°(B), which implies 

^ Tpdx — -^— f Tpdx — > \ilTp{z\) — p,lTp{z<2)- (4.19) 



Collecting l|4.16|l - H4.19|l . we find for it the equation 

Vu • Wlp dx + i/3 I uTpdx = fil (jp(z±) — Tp{z2)). 

This implies that u satisfies the first equation of Problem l|3.2(l on B. Thanks to Lemma 
13.11 the whole sequence (u e ) converges to u weakly in W 1,P (B) and strongly in L q (B), for 

l < q < 2p/{2- p). □ 

Let us conclude by some remarks: 

1. It is clear that the analysis carried out in this paper can be easily extended to the case 
where the physical properties /i and a are not constant. We shall however assume, in 
this case, that the magnetic permeability is a W 1 ' 00 function. This is necessary for H 2 
regularity results. 

2. The obtained results are generalizable to an arbitrary number of ("thick" or "thin") 
conductors. 

3. In the particular case where no "thick" conductor is present [i.e. Q,q = 0), the limit 
problem becomes 



-Au = nI(5 Zl -S Z2 ) 
Clearly, the solution of this equation is given by 



m 



1.2 



u{x) = !f logj^-J*}, xeR 2 . 

Z7T \X — Z\ | 



14 



References 

[1] H. Brezis, Analyse fonctionnelle, Masson, Paris, 1983. 

[2] L. Boccardo, T. Gallouet, Non-linear elliptic and parabolic equations involving measure 
data, J. Funct. Anal., Vol. 87 (1989) 149-169. 

[3] A. Damlamian, Ta Tsicn Li, Comportement limite des solutions de certains problemes 
mixtes pour les equations paraboliques, J. Math. Pures et Appliquees, 61 (1982) 113-130. 

[4] A. Friedman, Partial Differential Equations, Holt, Rinchart and Winston, Inc., 1969. 

[5] D. Gilbarg, N.S. Trudinger, Elliptic Partial Differential Equations of Second Order, 
Springer- Verlag, Berlin - Heidelberg - New York (1977). 

[6] J.-L. Lions, E. Magcncs, Problemes aux limites non homogenes et applications, Tome I, 
Dunod, Paris, 1968. 

[7] F. Murat, Soluciones renormalizadas de EDP elipticas non lineales, Preprints of the 
Laboratoirc Jacques-Louis Lions No. 93023, Universite Pierre et Marie Curie, 1993. 

[8] J. Necas, Les methodes directes en theorie des equations elliptiques, Masson, Paris (1967). 

[9] J.-C. Nedelec, Notions sur les equations integrates de la physique, Centre de 
Mathematiques Appliquees, Ecole Polytechnique, Palaiseau (1977). 

[10] J. Rappaz, R. Touzani, On a two-dimensional Magnetohydrodynamic problem, I. Mod- 
elling and analysis, Modelisation Mathematique et Analyse Numerique, Vol. 26, No. 2 
(1992) 347-364. 



15 



